Amplitude Modulation and Relaxation-Oscillation of Counterpropagating Rolls within 
a Broken-Symmetry Laser-Induced Electroconvection Strip 
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We report a liquid-crystal pattern-formation experiment in which we break the lateral (transla- 
tional) symmetry of a nematic medium with a laser- induced thermal gradient. The work is motivated 
by an improved measurement (reported here) of the temperature dependence of the electroconvec- 
tion threshold voltage in planar-nematic 4-methoxybenzylidene-4-butylaniline (MBBA) . In contrast 
with other broken-symmetry-pattern studies that report a uniform drift, we observe a strip of 
counterpropagating rolls that collide at a sink point, and a strong temporally periodic amplitude 
modulation within a width of 3-4 rolls about the sink point. The time dependence of the ampli- 
tude at a fixed position is periodic but displays a nonsinusoidal relaxation-oscillation profile. After 
reporting experimental results based on spacetime contours and wavenumber profiles, along with a 
measurement of the change in the drift frequency with applied voltage at a fixed control parameter, 
we propose some potential guidelines for a theoretical model based on saddle-point solutions for 
Eckhaus-unstable states and coupled complex Ginzburg-Landau equations. Published in PRE 73, 
036317 (2006). 



I. INTRODUCTION 

The formation of spatially periodic patterns in a wide 
variety of media driven far from equilibrium continues to 
attract a great deal of attention from experiment, simu- 
lations, and theory 0,|llllillll0,lllilI3IIlllllll 

For patterns in fluids, many important advances have 
relied on careful experiments with samples possessing a 
very high degree of spatial homogeneity in the material 
properties and the external control parameters. This 
homogeneity has proved crucial for careful comparisons 
between experiments and theory, since the latter usu- 
ally assumes that the medium has completely uniform 
spatial properties, resulting in a high degree of symme- 
try [lj|. Experimental patterns within samples of high 
purity and uniformity will undoubtedly continue to be 
of major importance in the foreseeable future. On the 
other hand, patterns in nature almost never occur in 
highly homogenized media [l5j , and one may reasonably 
inquire whether an inhomogeneous fluid, with a symme- 
try that is broken due to the variation of either a ma- 
terial property or a control parameter (or both) along 
some specific direction, will give rise to significant new 
physics. In the case of convection in horizontal fluid cells, 
broken up-down (U/D) symmetry about the midplane 
(non-Boussinesq convection) has been studied intensively 
for several decades and now forms a cornerstone of the 
Raylcigh-Bcnard convection (RBC) problem [llllj- Al- 
though broken lateral symmetries (i.e., broken transla- 
tional symmetry) in fluid media have received somewhat 
less attention than broken U/D symmetry, important re- 
sults have been obtained. A seminal work for broken lat- 
eral symmetry was the 1982 theoretical study by Kramer 
and coworkers |l6| on very small material-property ramps 
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in pattern-forming media, which motivated a number of 
new experiments and simulations on ramped systems. 
Hartung et al. \li\ . for example, discovered convection 
patterns with a uniform group drift in a RBC experiment 
in which both U/D and lateral symmetry were broken 
due to the temperature-dependence of the viscosity and 
spatial modulations in the cell height. Rehberg et al. p!q 
showed that two ramps with opposite slope could force 
uniform phase drift in RBC. Ning et al. [13 found travel- 
ing vortex waves in a Taylor vortex experiment in which 
a ramp was created with an axial variation of the annular 
width. Wiener and coworkers |2Cll | investigated ramped 
Taylor-vortex flow using an hourglass cell, which forms 
a particularly dramatic and intriguing example of bro- 
ken lateral symmetry. Beyond fluids, there has also been 
considerable interest in spatially varying parameters in 
reaction-diffusion Turing systems |2lj . 

It is the purpose of the present article to extend the 
understanding of pattern-formation in fluids with a bro- 
ken lateral symmetry. We report an experiment in which 
a controlled spatial ramp of the temperature within a 
nematic liquid crystal (NLC) produces spatially localized 
electroconvection rolls subjected to a positive nonuniform 
control parameter. In a planar-aligned NLC, a preferred 
direction for orientation of the pattern is selected by the 
nematic director. We create a long narrow strip-pattern 
parallel to the director, and break the lateral (transla- 
tional) symmetry along the strip using a laser-driven hor- 
izontal thermal gradient. The liquid crystal is planar- 
nematic 4-methoxybcnzylidene-4-butylaniline (MBBA). 
Our experiments, which expand on a preliminary study 
limited to much smaller pattern areas |22], rely on (1) 
the use of nonuniform laser irradiation to produce a long 
strip parallel to the director that is slightly warmer than 
its surroundings, and (2) the fact that the control pa- 
rameter in this medium is (in general) strongly depen- 
dent on temperature, as we show directly in experimen- 
tal detail below. The experimental results are of interest 



because they differ qualitatively both from known near- 
threshold electroconvection patterns in planar MBBA 
[23,121,123 J ^^^ from work (cited above) on other systems 
with a broken lateral symmetry. We find that the tem- 
poral variation of the roll amplitude at a fixed position is 
clearly nonsinusoidal. We investigate the new dynamics 
by studying spacetime diagrams, the temporal variation 
of the amplitude at fixed positions, and wavenumber pro- 
files. Although a theory applicable to electroconvection 
in liquid crystals with strong spatial inhomogeneities is 
not available at present, we proffer suggestions, based 
on the experimental results, for potentially useful ap- 
proaches to developing a theoretical model to explain 
our observations. Many of these considerations should 
be applicable to other broken-symmetry pattern-forming 
fiuids driven outside of equilibrium. 



II. RESULTS AND DISCUSSION 

The experimental apparatus and samples employed for 
our measurements have been described previously |22l- 
Briefly, a nematic film of MBBA (Aldrich) is doped with 
a small amount of blue-absorbing dye (methyl red, MR) 
and placed within a conventional homemade electrocon- 
vection cell with a thickness of 40 /ini and illuminated 
from above with a 488-nm laser beam. The electrocon- 
vection pattern is viewed using an inverted microscope 
and a CCD camera with a 640 X 480 pixel array and a 
resolution of 0.54 pixels//^m. We choose spatial coordi- 
nates such that the x and y axes in the nematic plane 
are parallel and perpendicular to the nematic director, 
respectively. For quantitative measurements of the pat- 
tern, we employ the background-divided dimensionless 
intensity, defined as |3ll29| 
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where G(x,y,t) and Gb{x,y) are the sample and back- 
ground gray-scale values, respectively, obtained at the 
pixel (x, y) of the CCD array at time t. The sample and 
background images are obtained under identical condi- 
tions except that for the latter, the voltage V across 
the sample is set to zero. At very low applied frequen- 
cies, we observed that the dimensionless intensity was 
spatially periodic along the director with a period equal 
to the sample thickness, as expected from shadowgraph 
theory [S^, ISll l33 |. The spatial period of the dimen- 
sionless intensity is thus one-half the spatial period of 
the director-tilt angle that defines the amplitude of the 
pattern. The control parameter for an electroconvection 
experiment is defined as e = V^ /Vq — 1, where Vq is the 
critical threshold voltage at which rolls first appear. In 
our experiments, e varies with position because of the 
non-uniform laser heating. We consistently measure the 
threshold voltage Vq at which rolls first appear along the 



laser-illuminated strip as the voltage is increased from 
zero. As explained below, the rolls occur first at the 
point at which the laser intensity is largest. Therefore, 
when the voltage across the sample is raised to y > Fq, 
the difference V'^/Vq — 1 represents the maximum emax 
in the control parameter as a function of position. We 
will employ emax defined in this way as our experimental 
control parameter. 

If the control parameter e that drives pattern formation 
in a fiuid is different at different temperatures, one may in 
principle create controlled spatial variation in the pattern 
structure using thermal inhomogeneities. We present in 
Fig.^measurements obtained at various applied frequen- 
cies showing the MBBA EC threshold voltage as a func- 
tion of temperature. No laser beam is present for these 
experiments, and the threshold voltage at a given applied 
frequency is simply defined as the voltage at which rolls 
are first observed on the sample at a fixed point that is 
far from the edges. A nonmonotonic temperature depen- 
dence was also reported in our previous study i2^] at a 
single applied frequency of 70 Hz; we see from the cur- 
rent results that, away from the minimum in Vo(T), the 
slope magnitude | dVo/dT \ increases sharply with applied 
frequency. These results are consistent with the very 
recent threshold-voltage measurements of Toth-Katona 
and Gleeson Q on the nematic Mischung V (M5), which 
show that I dVo/dT \ increases with frequency. In the 
M5 data reported by Toth-Katona and Gleeson, dVo/dT 
is always negative. Toth-Katona and Gleeson carefully 
noted that their EC patterns had greater spatial inho- 
mogeneity at higher frequencies. In the present work, 
as with our earlier study, we exploit this type of inhomo- 
geneity to control the spatial variation of the pattern am- 
plitude with an applied symmetry-breaking lateral tem- 
perature gradient. Indeed, from the data of Fig. ^ one 
may envision two methods for controlled variation in pat- 
tern structure using a region that is warmer than its sur- 
roundings. (1) At a bulk sample temperature such that 
dVa/dT < 0, one can generate a locally warmed region 
with e > surrounded by a linearly stable conducting 
background with e < 0, resulting in a spatially local- 
ized pattern. (2) If instead dVo/dT > at the bulk 
temperature, then one can warm a local region within a 
preexisting pattern to generate a local amplitude "hole" 
which (relative to the surrounding pattern) has a smaller 
e and hence a reduced amplitude. As discussed below, 
we have used laser light to successfully generate a local 
region with e > and hence a localized pattern (pro- 
cess (1)). We have found, however, that process (2) is 
much more difficult experimentally, because according 
to Fig. n] the required bulk temperature is close to the 
nematic-isotropic transition temperature, and we often 
observed the formation of small regions of isotropic fluid 
that strongly distorted the pattern. Hence in this article 
we will confine our attention to localized patterns pro- 
duced via process (1) {i.e. dVo/dT < 0). In passing, 
we should also note that if one wishes to minimize the 
effects of small thermal inhomogeneities in MBBA, one 
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FIG. 1: Measurement of the electroconvection threshold voh- 
age in MBBA as a function of temperature, relative to the 
nematic-isotropic transition temperature Tni, carried out at 
different applied frequencies. Each curve is labeled with the 
applied frequency in Hz. Data points are connected with lines 
to guide the eye. 



should work at the temperature at which Vo{T) achieves 
a minimum at the operating frequency. 

Our experimentally observed nonmonotonic tempera- 
ture dependence for the threshold voltage in MBBA can 
be compared to the theoretical "Standard Model" for EC 
develo ped by Bodenschatz, Kramer, Pesch, and Zimmer- 
mann |33 . l34Ll35l |. In this model, the square of the thresh- 
old voltage Vq can be expressed as stated by Kramer and 
Pesch HJI: 
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where K^^^f\ea ,(7a , and rj^^ff^ represent effective 
parameters for the elasticity, the dielectric anisotropy, 
the anisotropy in the electrical conductivity, and the 
viscosity, respectively. For nonobliquc rolls, each of 
these four effective parameters depends on the wavenum- 
ber q selected by the pattern; furthermore, ea and 
aa also depend on the frequency of the applied 

field. The constants a2 and Tq represent the second 
Leslie coefficient and the charge relaxation time, respec- 
tively. The temperature dependence of each parameter 
in E g. [3 can be evaluated using studies from the litera- 
ture [Sa, I33, lla, I33 . We normalized the literature val- 
ues of the electrical conductivity to our measured room- 
temperature MBBA value of cr^ = 8 x 10-^(f}m)-i [2l |. 
At each temperature T, the correct wavenumber q was 
determined with a spreadsheet based on Eq.[21by adjust- 



FIG. 2; Prediction of the threshold voltage at different tem- 
peratures using the Standard Model. The open and filled cir- 
cles correspond to applied frequencies of 70 Hz and 200 Hz, 
respectively. 



ing q until Vq achieved a minimum, which provided the 
desired Vq{T). As an independent test of our spread- 
sheet algorithm, we calculated Vq using as inputs the 
room-temperature MBBA parameters published by Bo- 
denschatz et al. [33. Figure 3 of that paper shows Vq 
vs. applied frequency, based on the Standard Model with 
rigid boundary conditions at 25 °C. We found that our 
spreadsheet-calculated values for both the selected q and 
the threshold voltage Vq at ten different applied frequen- 
cies agreed with those of Bodenschatz et al. at a level 2% 
or better. Our spreadsheet-based Standard-Model results 
for Vo as a function of temperature are shown in Fig. |21 
at two different applied frequencies. The nonmonotonic 
temperature dependence and the strong increase in the 
slope I dVo/dT \ with applied frequency are both apparent 
in the theoretical Standard-Model prediction depicted in 
Fig. [3 The values for Vq predicted using the theory do 
not agree quantitatively with the experiment: although 
the theory and the experimental data are consistent to 
about 20% at 70 Hz, the difference is about a factor of 
two at 200 Hz. This is not a serious concern, since there 
are four separate effective material parameters in the the- 
ory, each of which (most notably the electrical conductiv- 
ity) displays considerable variation in the literature. The 
tendency of Vo{T) to increase at both the high and low 
temperature limits of the nematic phase is quite clear in 
the Standard Model, and we believe our experiments are 
consistent with this model. 

The temperature dependence of the threshold voltage 
in MBBA reported in Fig.^above confirms that it should 
be possible to use lateral temperature gradients in MBBA 
to create interesting variations in EC pattern structure, 
which is the basis of the experiments reported below. 



These measurements are an improvement on a previous 
study |23| in which eUiptical locaHzed EC pulses with 
major axes of 30 - 300 /im were created in MR-doped 
MBBA using a laser beam with a circular cross-section 
with a diameter of about 50 /im. In that work we con- 
cluded that the mechanism for the pulse creation was 
thermal, because essentially identical patterns were seen 
using anthraquinone as an alternative absorbing dye. In 
the previous work, we found that the amplitude and lat- 
eral size of the laser-induced convection pattern actu- 
ally decreased when the control parameter was driven 
past emax ^ 0.4. This unexpected attenuation was at- 
tributed to the advection-driven increase in the effec- 
tive thermal conductivity (or equivalently, an increase in 
the local Nusselt number) at larger pattern amplitudes, 
which quenches the temperature spike produced by the 
absorbed laser light. We also found interesting roll drift 
in regions where the laser intensity was spatially nonuni- 
form, but not within regions in which the laser intensity 
was uniform. In the current study, we have made two 
experimental improvements, resulting in a simplification 
in the pattern structure and the dynamics. First, we 
have used cylindrical lenses to stretch the laser cross- 
section into a long strip parallel to the director. For a 
bulk sample temperature at which dVo/dT < 0, a laser- 
light strip of this type can be used to produce a quasi- 
one-dimensional pattern delimited by "soft" |l3l | bound- 
ary conditions: a long narrow strip in which e{x) > is 
surrounded by a linearly stable state with e < 0. Long 
narrow electroconvection channels have been generated 
previously using etching techniques and "hard" bound- 
ary conditions by Hidaka et al. |40j, and Rasenat, Braun, 
and Steinberg |4lLl43|. Second, we have fixed the control 
parameter at emax =0.1 throughout the current study, to 
avoid complications caused by the amplitude dependence 
of the local Nusselt number. (In an EC experiment, the 
constraint Smax = 0.1 maps out a curve that is displaced 
slightly from the usual emax = threshold curve in the 
plane of applied frequency vs. applied voltage.) To pro- 
duce a nonuniform laser-intensity profile, we used a razor 
blade to block one-half of the length of the rectangular 
cross-section, so that knife-edge diffraction resulted in 
an inhomogeneous laser-beam profile along the director. 
A typical laser intensity profile is shown in Fig. |2| By 
measuring the change in the nematic-isotropic transition 
temperature at the point of maximum intensity within 
the laser-illuminated region, using both increasing and 
decreasing bulk temperatures, we inferred that the tem- 
perature increase at this point was about 0.9 °C. This 
temperature increase is not large enough to significantly 
change the dye absorption coefficient. 

Images of a localized pattern of short rolls produced 
along the laser strip at three different times with a control 
parameter of €max = 0.1 are shown in Fig.^ In contrast 
to previous patterns commonly observed in MBBA close 
to threshold, in which the rolls are either stationary or 
drift slowly in one direction |2a.l2l.l25l |. the rolls localized 
along the laser strip show phase drifts that counterprop- 
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FIG. 3; A typical laser intensity profile. The horizontal and 
vertical axes represent x and y, respectively, in units of mi- 
crometers. The width in the y direction (FWHM) is about 
11 ^m. The total laser power incident on the sample for 
the experiments was 70 ^W. A HIGHER-RESOLUTION 
PLOT OF THE LASER PROFILE IS AVAILABLE IN THE 
PUBLISHED VERSION OF THIS DOCUMENT (PRE 73, 
036317, 2006). 






FIG. 4: Direct images (without background correction) of the 
electroconvection patterns formed along the laser-light strip 
at several different times. From bottom to top, the images 
were acquired at times f = 0, 9 s, and 27 s. The applied 
voltage and frequency were 9.33 V and 200 Hz, respectively, 
which resulted in a control parameter of emax =0.1. The bulk 
sample temperature was T — Tni = —16.7 °C. The white bar 
in the lower-left corner represents a length of 100 fira, and 
the short white vertical bar to the right of center in the lower 
panel indicates the position x = 0. The image contrast has 
been enhanced so that the low-amplitude rolls are easier to 
see. 



agate into a sink point {x = 0), as shown with spacetime 
contours in Fig. [3 As the rolls drift, the ft'equency froii of 
the amplitude variation at a specific position x was found 
to be constant with respect to x to within about 1%. In 
Fig. spacetime contours for the same images at higher 
spatial and temporal resolution show that in addition to 
counterpropagating phase drifts, there is also a strong 
temporally periodic amplitude modulation. Indeed, an 
amplitude "gap" in a central region about the sink point 
periodically appears in which 3-4 rolls are reduced to es- 
sentially zero amplitude. (We will subsequently refer to 
the region in which the amplitude gap forms as the "cen- 
tral zone".) A new set of rolls then rapidly grows into 
the gap region, until the amplitude near .t = is actu- 
ally higher than the amplitude outside the central zone. 
Next, the central-zone rolls slowly decrease in amplitude 
until a gap has again formed. This new gap is then again 
replaced by a set of rolls, but the new rolls are shifted 
in phase by about one-half cycle relative to the rolls that 
grew into the previous gap. For the sample employed for 
Fig.|^ we varied the applied frequency over the range / = 
160 — 430 Hz, while simultaneously varying the applied 
voltage over the range V = 8.2 — 57 V to keep the con- 
trol parameter fixed at emax = 0.1. We examined space- 
time contours at fourteen different settings for (/, V) and 
found that the same basic dynamics, with counterprop- 
agating phase drifts and a strong amplitude modula- 
tion leading to a central-zone gap, was clearly present 
at each (/, V) setting. The speed of both the phase 
drift and the amplitude modulation was found to increase 
strongly with applied voltage, as we report below in de- 
tail. In addition, this same basic counterpropagating- 
phase-drift/amplitude-modulation dynamic has been ob- 
served on four different independently prepared samples 
illuminated with a laser intensity that was nonuniform 
in the x-direction. The dashed line in Fig. |51 shows 
that in addition to the counterpropagating phase drift, 
there is also a slow overall group drift of all the rolls to- 
ward positive X. From the figure one may conclude that 
the group velocity for overall drift is about 1/7 of the 
phase velocity corresponding to the counterpropagating 
roll drift. An overall group drift is commonly observed 
in the literature for patterns subjected to a broken lat- 
eral symmetry 0, ^ |^ . On the other hand, the 
counterpropagating-roU phase drift and the temporally 
periodic amplitude modulation in the central zone are 
unexpected, and will require further elucidation. 

From Fig.Elit is clear that the temporal dependence of 
the amplitude at a fixed point x within the central zone 
is periodic but not sinusoidal. This is shown more explic- 
itly in Fig. [3 which displays the amplitude at x = for 
the data of Fig. 0] as a function of time for one period at 
two different applied voltages. The rising edge is consid- 
erably faster than the falling edge. This type of asym- 
metric temporal behavior, which is usually referred to 
as relaxation-oscillation in the nonlinear dynamics liter- 
ature |43| . is fundamentally nonlinear, since linear insta- 
bility analysis will always result in a (in general complex) 
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FIG. 5: A spacetime plot showing the evolution of the pat- 
tern displayed in Fig. |31 Each trace represents, at different 
times, the dimensionless intensity averaged over the y direc- 
tion. Time evolves upward, and the total elapsed time is 
450 s. 



exponential time dependence. 

As the phase-drift and amplitude-modulation pro- 
cesses occur, the wavenumber profile q{x) of the pattern 
varies with time. The variation is shown in Fig.|Sl where 
q{x), normalized to its value for the left-most spatial pe- 
riod in the pattern, is shown at several different times 
before the amplitude gap forms. At a given time, the 
wavenumber q(x) has its maximum near x ~ 0, and is 
nearly symmetric about that maximum. The maximum 
fractional difference between the wavenumber near the 
left edge of the strip and the wavenumber at x = is 
about 30%. This exceeds the Eckhaus bandwidth for 
stable periodic states at a control parameter of e — 0.1, 
which is \/e/3 ~ 18% |4j. We may therefore antici- 
pate that wavelength-changing processes driven by the 
Eckhaus instability are potentially relevant in our exper- 
iments. Our observation in Fig. that the location of 
the amplitude the gap is selected by the system in a cen- 
tral region about a; = (i.e., far from the boundaries) 
appears to be consistent with predictions of Kramer and 
Zimmermann |2y| regarding the location of the reduced- 
amplitude wavelength-changing zone within an Eckhaus- 
destabilized pattern. These authors argued that, while 
the location of the wavelength-changing process in an in- 
finite system is arbitrary, in a real system this process 
will occur in a centralized "weak" region far from the 
boundaries, in a manner consistent with Tay lor vortex 
flow experiments |43 and simulations [lj,|4g. 

It is also of interest to compare the roll-dynamic 
timescales that we observe experimentally to the relevant 
Giznburg-Landau (GL) timescale. The one-dimensional 
GL equation for a slowly-varying amplitude A{x, t) in the 
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FIG. 6: A spacetime plot of the pattern at higher spatial and 
temporal resolution. The region referred to as the "central 
zone" in the text corresponds to the width of the amplitude 
gap, which is indicated (at a specific time) with the two plus 
signs. Six traces showing the pattern that grows into the gap 
are indicated with thick lines; note that alternate thick-line 
traces are phase-shifted by about one-half of a cycle. The 
slow overall group drift is indicated by the dashed line, which 
has a reciprocal slope equal to the group velocity. 



FIG. 8: Variation of the wavenumbcr with position along the 
X-axis. The data employed are the earliest six traces from the 
spacetime plot of Fig. |S1 before the first amplitude gap forms. 
The wavenumber is normalized to its value on the far left. 
Time evolves upward, with an interval of 3 s between consec- 
utive plots. The thick vertical bar on the lower part of the 
left axis represents a wavenumber variation of 30%; thus, we 
observe on this figure that the peak normalized wavenumber 
near x = increases (with time) from about 1.1 to about 1.3. 



weakly nonlinear regime (e << 1) can be written [T3lll4 | 
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FIG. 7: Measured dimensionless intensity (see Eq. as a 
function of time at the position a:: = 0, using two differ- 
ent applied voltages/frequencies. The control parameter is 
Smax = 0.1 for both curves. The applied frequency and volt- 
age is (400 Hz, 31.4 V) for the lower curve, and (200 Hz, 
9.33 V) for the upper curve. The total elapsed time along the 
horizontal axis is 4.3 s (lower) and 36 s (upper). The upper 
curve has been offset vertically by one unit. 



TdtA{x,t) = ^;diA{x,t)+eA{x,t)-go \ A 
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in which the space and time scales are represented by 
^0 and T, respectively, with e serving as the control 
parameter and go as the usual negative-feedback coeffi- 
cient. For clectroconvcction, r is given approximately by 
the director relaxation time Tdir, since this is generally 
the longest timescale appearing in the microscopic 
equations of motion |33l |3J|. It is possible that the 
control parameter e can vary on position and time scales 
that are large compared to ^o and r, respectively |47| . 
To compare the GL timescale to our experimentally 
measured timescales, we designed an experiment in 
which we measured the temporal frequency froii of the 
intensity variation at a fixed point x as a function of the 
applied voltage across the sample. In this measurement, 
it is again important to note that the applied voltage and 
frequency were varied together in a manner such that 
the control parameter was held constant at emax = 0.1. 
Results for the measured frequency froii vs. the applied 
voltage across the sample are shown in Fig. |51 with frou 
measured at x = 0. By using the Standard Model to 
compute Tdir at each applied voltage, we may also show 
^max/Tdir, scalcd down by a factor of 10~^, on the same 
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FIG. 9: The individual points are the measured roU frequen- 
cies froii at the position a; = shown as a function of the 
square of the apphed voltage V. The applied frequency and 
voltage were both varied so that the control parameter e,nax 
is equal to 0.1 at each individual point displayed. The solid 
line is the director relaxation rate e,nax/Tdir, scaled by a fac- 
tor of lO"'^ so that froii and emax/Tdir can be displayed on 
the same graph. 



group drift results from the thermal ramp. However, the 
observed counterpropagating phase drift and the central- 
zone amplitude modulation do not appear to be compati- 
ble with previously reported drift mechanisms in ramped 
media for two reasons. First, although our ramp, in the 
form of an imposed thermal gradient, has the same sign 
everywhere in the pattern, the observed phase drift is 
counterpropagating rather than unidirectional. Second, 
in addition to phase drift, we observe a strong periodic 
amplitude modulation in the region close to x = 0. The 
incompatibility of our results with ramp-theory drifts is 
perhaps not surprising in view of the requirement of very 
small ramps in the theory. 

To the best of our knowledge, all theoretical work on 
models of electroconvection has been carried out assum- 
ing a uniform isothermal medium for the spatially ex- 
tended rolls. Thus, a theoretical model of electroconvec- 
tion patterns subjected to a local temperature gradient, 
and in particular the use of local heating to produce lo- 
calized pattern formation with soft boundary conditions, 
has not yet been formulated. Based on the results above, 
we believe that, in closing, we can offer suggestions on 
potentially useful avenues for theoretical development on 
these interesting dynamics. In developing a new model, 
one must demand explanations for the following salient 
aspects of the data, which represent distinct new fea- 
tures vis-a-vis traditional spatially extended EC rolls in 
uniform MBBA: 



graph. Clearly, at a given applied voltage, the measured 
frequency froii is always much smaller than tmax/Tdir- 
Thus the temporal variations in amplitude and phase 
that we observe is a slowly varying temporal modulation 
vis-a-vis the GL equation. 



• opposing phase drifts for very negative and very 
positive X, resulting in counterpropagating rolls; 

• a strong amplitude modulation resulting in the 
temporally periodic formation of a "gap" in the 
central region, with a temporal profile that is peri- 
odic but not sinusoidal. 



III. CONCLUSIONS AND FUTURE WORK 

In considering potential theoretical models to explain 
our experiments, wc consider first the theory of very 
small parameter ramps in pattern- forming media due 
to Kramer and coworkers [l6j. A ramp within a one- 
dimensional pattern selects a specific (though in general 
nonuniversal) wavenumber profile function q{x). It is well 
known that drift can occur in ramped medium. Riecke 
and Paap [43, for example, showed that if q{x) pene- 
trates into the Eckhaus-unstable baiid, phase drift may 
occur. In addition, Rchberg et al. |18J carried out ex- 
periments in which two ramps with opposite slope at 
the ends of a Rayleigh-Bernard cell selected two differ- 
ent wavenumber profiles qieft{x) and qright{x), which can 
generate a unidirectional phase drift from high q toward 
low q. The authors explained these experiments with 
a qualitative model showing that drift will occur when- 
ever there is no overlap between the finite- width selection 
bands centered on qieft{x) and qright{x). For our experi- 
ments, it is reasonable to conclude that the slow uniform 



We wish to point out two theoretical approaches that 
could prove useful in developing a model for these exper- 
imental features. 

(1) First, it is clear from Fig. |H1 that the wavenumber 
bandwidth is driven back and forth across the Eckhaus 
boundary. Therefore, if one wishes to apply solutions 
to the GL equation applicable to our experiments, one 
should include not only the usual constant-amplitude so- 
lutions A{x, t) = gg y^e — k^^Q e^^^ , where k = q — qo, 
but also the ampfitude-modulated saddle-point (SP) so- 
lutions that separate basins of attraction of the constant- 
amplitude solutions. These SP solutions, which drive the 
wavelength-changing processes that arise when the Eck- 
haus boundary is crossed, were considered in detail in 
two widely cited papers by Kramer, Schobcr, and Zim- 
merman [2^ I23 . The simplest SP solution considered by 
these authors is the function 



Asp{x,t) 



Akx 



y/2 k^o + iS{t) tanh( 






(4) 



where S{t) = \Je{t) — Sfc^^Q is real, and we have in- 
cluded a potentially slowly- varying control parameter 
e(i). Kramer and Zimmermann [26j pointed out that 
for large | a; | , Asp takes the form of the usual periodic 
GL solution (proportional to e**^^); however, due to the 
tanh function, there is a phase shift between the periodic 
solutions at x = ± oo. In view of the experiments, an 
interesting characteristic of Eq. 0] is the fact that if the 
control parameter is changed slowly, the phase of Asp 
will change in opposite directions at a; = ± oo. 

2. Since counterpropagating traveling waves are seen 
under all experimental conditions explored, it is certainly 
reasonable to consider application of two coupled com- 
plex Ginzburg-Landau equations (GGLE)[43|, in which 
the parameters ^q and g^ in Eq. |21 above will have 
both real and imaginary parts. A direct comparison 
of experiment with CGLE-based theory can be difficult, 
because the CGLE can describe an extraordinary vari- 
ety of different traveling- wave dynamics, and the coef- 
ficients of the equation are in general difficult to deter- 
mine expcrimentallv[5(]j . However, it is well-known that 
a reduction of the dimensionality or localization in the 
experimental patterns, which have both been achieved 
in the present experiments by using a long strip of laser 
light to localize the rolls, can reduce considerably the 
difficulty of determining the coefficients. A useful start- 
ing point would appear to be the CGLE modeling ap- 
proach adopted by Pastur et a/.|50l| in describing their 
heated-wire experiments. These authors observed left- 
and right-traveling waves separated by sources and sinks. 
If the CGLE model employed by Pastur et al. could be 
extended to the case of nonuniform spatial heating by in- 
troducing a spatially dependent control parameter e(a;), 
it would be very interesting to see if the traveling-wave 
dynamics resulting from the model displayed character- 
istics that reproduced the experimental results above. 

To conclude, we have presented experimental results on 
electroconvection-pattcrn formation in a narrow strip in 
which the lateral symmetry has been broken with a laser- 
induced thermal gradient. Our main purpose was to re- 
port new experimental results that differ in a qualitative 
manner from previous experiments on broken-symmetry 
fluidizcd patterns. Our measured temperature depen- 
dence of the threshold voltage showed that is possible 
to create a thermally localized supercritical strip. The 
spacetime profiles along the strip displayed a slow coun- 



terpropagating phase drift and an asymmetric relaxation- 
oscillation temporal profile, in addition to a strong slowly 
varying temporally periodic amplitude modulation that 
generates a central-zone gap. Measurement of the spa- 
tial variation of the wavenumber displayed a wavenum- 
ber bandwidth that exceeded the Eckhaus stability limit. 
We noted some simple connections to Eckhaus unstable 
states and to the complex Ginzburg-Landau equation in 
an attempt to provide some potential guidance for the 
development of a theoretical model to explain the results. 
We are at present upgrading our apparatus and soft- 
ware to gain an improved understanding of the patterns 
discussed above. There are two future tasks that we be- 
lieve are especially exigent. First, it is very important 
to carry out a direct measurement of the temperature 
as a function of position along the pattern strip. This 
would allow a determination of the control parameter 
e(x) as a function of position, which is important for 
the development of a model. This is difficult experimen- 
tally since the maximum temperature variation is only 
on the order of 1 °C, but is crucial nonetheless. Second, 
one should investigate whether numerical simulations of 
coupled CGL equations with a driving force consistent 
with Fig.|21can show an amplitude modulation and coun- 
terpropagating phase drifts akin to the experimental re- 
sults. Code-writing for such simulations is currently in 
progress. 
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